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Abstract. Given a point A in the real Grassmannian, it is well-known that one can construct a 
soliton solution ua(x, y, i) to the KP equation. The contour plot of such a solution provides a tropical 
approximation to the solution when the variables x, y, and t are considered on a large scale and the 
time t is fixed. In this paper we give an overview of our work on the combinatorics of such contour 
plots. Using the positroid stratification and the Deodhar decomposition of the Grassmannian (and in 
particular the combinatorics of Go-diagrams), we completely describe the asymptotics of these contour 
plots when y or t go to ±00. Other highlights include: a surprising connection with total positivity and 
cluster algebras; results on the inverse problem; and the characterization of regular soliton solutions 
- that is, a soliton solution ua(x, y, t) is regular for all times t if and only if A comes from the totally 
non-negative part of the Grassmannian. 
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1. Introduction 

The main purpose of this paper is to give an exposition of our recent work [18, 19, 20], which found 
surprising connections between soliton solutions of the KP equation and the combinatorics of the real 
Grassmannian. The KP equation is a two-dimensional nonlinear dispersive wave equation which was 
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proposed by Kadomtsev and Peviashvili in 1970 to study the stability problem of the soliton solution 
of the Korteweg-de Vries (KdV) equation [13]. The equation has a rich mathematical structure, and is 
now considered to be the prototype of an integrable nonlinear dispersive wave equation with two spatial 
dimensions (see for example [24, 1, 9, 23, 12]). The KP equation can also be used to describe shallow 
water wave phenomena, including resonant interactions. 

An important breakthrough in the KP theory was made by Sato [26], who realized that solutions 
of the KP equation could be written in terms of points on an infinite-dimensional Grassmannian. The 
present paper, which gives an overview of most of our results of [19] and [20], deals with a real, finite- 
dimensional version of the Sato theory. In particular, we are interested in soliton solutions, that is, 
solutions that are localized along certain rays in the xy plane called line-solitons. Such a solution can 
be constructed from a point A of the real Grassmannian. More specifically, one can apply the Wronskian 
form [26, 27, 11, 12] to A to produce a certain sum of exponentials called a T-function TA{x,y,t), and 
from the T-function one can construct a solution UA{x,y, t) to the KP equation. 

Recently several authors have studied the soliton solutions ua(x, y, t) which come from points A of the 
totally non-negative part of the Grassmannian (Gr/ Si „)>o, that is, those points of the real Grassmannian 
Grk.n whose Pliicker coordinates are all non-negative [3, 16, 2, 5, 7, 18, 19]. These solutions are regular, 
and include a large variety of soliton solutions which were previously overlooked by those using the 
Hirota method of a perturbation expansion [12]. 

A main goal of [20] was to understand the soliton solutions ua(x, y, t) coming from arbitrary points A 
of the real Grassmannian, not just the totally non-negative part. In general such solutions are no longer 
regular - they may have singularities along rays in the xy plane - but it is possible, nevertheless, to 
understand a great deal about the asymptotics of such solutions, when the absolute value of the spatial 
variable y goes to infinity, and also when the absolute value of the time variable t goes to infinity. 

Two related decompositions of the real Grassmannian are useful for understanding the asymptotics 
of soliton solutions ua(x, y, t). The first is Postnikov's positroid stratification of the Grassmannian [25], 
whose strata are indexed by various combinatorial objects including decorated permutations and J- 
diagrams. This decomposition determines the asympototics of soliton solutions when |y| >• 0, and our 
results here extend work of [2, 5, 7, 18, 19] from the setting of the non-negative part of the Grassmannian 
to the entire real Grassmannian. The second decomposition, which we call the Deodhar decomposition 
[20], is the projection to the Grassmannian of a decomposition of the flag variety due to Deodhar. 
The Deodhar decomposition refines the positroid stratification, and its components may be indexed by 
certain tableaux filled with black and white stones called Go-diagrams, which generalize J-diagrams. 
This decomposition determines the asymptotics of soliton solutions when \t\ ^> 0. More specifically, it 
allows us to compute the contour plots at \t\ » of such solitons, which are tropical approximations 
to the solution when x, y, and t are on a large scale [20]. 

By using our results on the asymptotics of soliton solutions when (CO, one may give a characteri- 
zation of the regular soliton solutions coming from the real Grassmannian. More specifically, a soliton 
solution ua{x, y, t) coming from a point A of the real Grassmannian is regular for all times t if and only 
if A is a point of the totally non-negative part of the Grassmannian [20] . 

The regularity theorem above provides an important motivation for studying soliton solutions coming 
from (Grfc,„)>o. Indeed, as we showed in [19], such soliton solutions have an even richer combinatorial 
structure than those coming from Gr^^. For example, (generic) contour plots coming from the totally 
positive part (Gr^„)>o of the Grassmannian give rise to clusters for the cluster algebra associated to 
the Grassmannian. And up to a combinatorial equivalence, the contour plots coming from (Gr2. n )>a 
are in bijection with triangulations of an n-gon. Finally, if either A e (Gr fe-Il ) >0 , or A e (Gr fc .„)> and 
t <C 0, then one may solve the inverse problem for UA(x,y,t): that is, given the contour plot C t {uA) 
and the time t, one may reconstruct the element A e {Grk, n )>o- And therefore one may reconstruct 
the entire evolution of this soliton solution over time [19]. 
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The structure of this paper is as follows. In Section 2 we provide background on the Grassmannian 
and some of its decompositions, including the positroid stratification. In Section 3 we describe the 
Deodhar decomposition of the complete flag variety and its projection to the Grassmannian, while 
in Section 4 we explain how to index Deodhar components in the Grassmannian by Go-diagrams. 
Subsequent sections provide applications of the previous results to soliton solutions of the KP equation. 
In Section 5 we explain how to produce a soliton solution to the KP equation from a point of the real 
Grassmannian, and then define the contour plot associated to a soliton solution at a fixed time t. In 
Section 6 we use the positroid stratification to describe the unbounded line-solitons in contour plots of 
soliton solutions at y » and y <C 0. In Section 7 we define the more combinatorial notions of soliton 
graph and generalized plabic graph. In Section 8 we use the Deodhar decomposition to describe contour 
plots of soliton solutions for t -C 0. In Section 9 we describe the significance of total positivity to soliton 
solutions, by discussing the regularity problem, as well as the connection to cluster algebras. Finally in 
Section 10, we give results on the inverse problem for soliton solutions coming from (Gr fc .„)> . 

2. Background on the Grassmannian and its totally non-negative part 

The real Grassmannian Gr k , n is the space of all fc-dimensional subspaces of R n . An element of Gr k . n 
can be viewed as a full- rank k x n matrix modulo left multiplication by nonsingular k x fc matrices. In 
other words, two fc x n matrices represent the same point in Grk, n if an d only if they can be obtained 
from each other by row operations. Let be the set of all fc-element subsets of [n] :— {1, . . . , n}. For 
/ £ ('fc')) let A/(^4) be the Pliicker coordinate, that is, the maximal minor of the k x n matrix A located 
in the column set /. The map A h-» (Aj(A)), where / ranges over ('^), induces the Pliicker embedding 

Grk.n w( fc ) _1 . The totally non-negative part of the Grassmannian (Gr k , n )>o is the subset of Gr k . n 
such that all Pliicker coordinates are non-negative. 

We now describe several useful decompositions of the Grassmannian: the matroid stratification, the 
Schubert decomposition, and the positroid stratification. When one restricts the positroid stratification 
to (Gr k . n )>o, one gets a cell decomposition of (Grk, n )>o into positroid cells. 

2.1. The matroid stratification of Gr^„. 

Definition 2.1. A matroid of rank k on the set [n] is a nonempty collection M C (^) of fc-element 

subsets in [n], called bases of M, that satisfies the exchange axiom: 

For any I,J £ M and i £ I there exists j £ J such that (7 \ {i}) U {j} £ M. 

Given an element A £ Gr k , n , there is an associated matroid Ma whose bases are the fc-subsets 
I c [n] such that A/ (A) ^ 0. 

Definition 2.2. Let M C (^') be a matroid. The matroid stratum Sm is defined to be 

S M ={A£ Gr k . n | A 7 (A) ^ if and only if I £ M}. 

This gives a stratification of Gr kiTL called the matroid stratification, or Gelfand-Serganova stratification. 
The matroids Ai with nonempty strata Sm are called realizable over R. 

2.2. The Schubert decomposition of Grfe ; „. Recall that the partitions A C (n — k) k are in bijection 
with fc-element subset I C [n] . The boundary of the Young diagram of such a partition A forms a lattice 
path from the upper-right corner to the lower-left corner of the rectangle (n — k) . Let us label the n 
steps in this path by the numbers 1, . . . , n, and define I = as the set of labels on the fc vertical 
steps in the path. Conversely, we let X(I) denote the partition corresponding to the subset /. 

Definition 2.3. For each partition A C (n — k) k , one can define the Schubert cell £l\ by 

Q\ = {A £ Gr k n | /(A) is the lexicographically minimal base of Ma}- 
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As A ranges over the partitions contained in (n — k) , this gives the Schubert decomposition of the 
Grassmannian Gr kn , i.e. 

Gr k , n = \_\ ^A- 

\<Z(n-k) k 

We now define the shifted linear order <j (for i G [n]) to be the total order on [n] defined by 
i <i i + 1 <i i + 2 <i ■ ■ ■ <i n <i 1 <i ■ ■ ■ <i i — 1. 
One can then define cyclically shifted Schubert cells as follows. 

Definition 2.4. For each partition A C (n — k) k and i G [n], the cyclically shifted Schubert cell £l\ is 

£l\ = {A € Grk.n | 7(A) is the lexicographically minimal base of Ma with respect to <j}. 

2.3. The positroid stratification of Gr k ^ n . The positroid stratification of the real Grassmannian 
Gr kn is obtained by taking the simultaneous refinement of the n Schubert decompositions with respect 
to the n shifted linear orders <j. This stratification was first considered by Postnikov [25], who showed 
that the strata are conveniently described in terms of Grassmann necklaces, as well as decorated per- 
mutations, (equivalence classes of) plabic graphs, and d-diagrams. Postnikov coined the terminology 
positroid because the intersection of the positroid stratification with (G , r / t j „)>o gives a cell decomposi- 
tion of (Grfc jT j)>o (whose cells are called positroid cells). 

Definition 2.5. [25, Definition 16.1] A Grassmann necklace is a sequence X = (I±, . . . ,I n ) of subsets 
I r C [n] such that, for i S [n], if i € Ii then = (7j \ {i}) U {j}, for some j G [n); and if i ^ I r then 
7j + i = Jj. (Here indices i are taken modulo n.) In particular, we have |7i| = ■ ■ ■ = |7„|, which is equal 
to some k G [n]. We then say that X is a Grassmann necklace of type (k,n). 

Example 2.6. (1257, 2357, 3457, 4567, 5678, 6789, 1789, 1289, 1259) is a Grassmann necklace of type (4, 9). 

Lemma 2.7. [25, Lemma 16.3] Given A G Gr k ^ n , let 1(A) = (7i, . . . ,I n ) be the sequence of subsets in 
[n] such that, for i G [n], h is the lexicographically minimal subset of (^') with respect to the shifted 
linear order <, such that Ai^A) ^ 0. Then 1(A) is a Grassmann necklace of type (k,n). 

If A is in the matroid stratum 5m, we also use to denote the sequence (7i, . . . , 7„) defined above. 
This leads to the following description of the positroid stratification of Grk, n - 

Definition 2.8. Let I = (7i, . . . , I n ) be a Grassmann necklace of type (k,n). The positroid stratum Si 
is defined to be 

Si = {A G Gr fc ,„ | J(A) = 1} 

n 

»=i 

The second equality follows from Definition 2.8 and Definition 2.4. Note that each positroid stratum 
is an intersection of n cyclically shifted Schubert cells. 

Definition 2.9. [25, Definition 13.3] A decorated permutation w = (n, col) is a permutation ir G S n 
together with a coloring function col from the set of fixed points {i \ ir(i) — i} to {1, —1}. So a decorated 
permutation is a permutation with fixed points colored in one of two colors. A weak excedance of w is 
a pair (i, 7Tj) such that either n(i) > i, or ir(i) = i and col(i) = 1. We call i the weak excedance position. 
If n(i) > i (respectively ir(i) < i) then (i,7Tj) is called an excedance (respectively, nonexcedance) . 

Example 2.10. The decorated permutation (written in one-line notation) (6,7,1,2,8,3,9,4,5) has no 
fixed points, and four weak excedances, in positions 1, 2, 5 and 7. 
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Definition 2.11. A plabic graph is a planar undirected graph G drawn inside a disk with n boundary 
vertices 1, ...,n placed in counterclockwise order around the boundary of the disk, such that each 
boundary vertex i is incident to a single edge. 1 Each internal vertex is colored black or white. See the 
left of Figure 1 for an example. 

Definition 2.12. [25, Definition 6.1] Fix k, n. If A is a partition, let Y\ denote its Young diagram. A 
d-diagram (A, D)k, n of type (k, n) is a partition A contained in a k x (n — k) rectangle together with a 
filling D : Y\ — > {0, +} which has the J-property: there is no which has a + above it and a + to its 
left. 2 (Here, "above" means above and in the same column, and "to its left" means to the left and in 
the same row.) See the right of Figure 1 for an example. 
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k = 6, n= 16 

K = (10, 9, 9, 8,5,2) 



Figure 1. A plabic graph and a Le-diagram L = (A, Z?)fe, n . 

We now review some of the bijections among these objects. 

Definition 2.13. [25, Section 16] Given a Grassmann necklace X, define a decorated permutation 7T = 
7r(I) by requiring that 

(1) if = (h \ {»}) U {j}, for j i, then 7T : (j) = i. 3 

(2) if 2j_|_i = Ii and i G Zj then 7r(i) = i is colored with col(i) = 1. 

(3) if = Ii and i ^ Ii then = i is colored with col(i) = — 1. 

As before, indices are taken modulo n. 

If 7T = 7r : (I), then we also use the notation S v -- to refer to the positroid stratum Sj. 

Example 2.14. Definition 2.13 carries the Grassmann necklace of Example 2.6 to the decorated permu- 
tation of Example 2.10. 

Lemma 2.15. [25, Lemma 16.2] The map X — ¥ 7r(Z) is a bijection from Grassmann necklaces X = 
[I\ ,...,/„) of size n to decorated permutations 7r (I) of size n. Under this bijection, the weak excedances 
of 7T (X) arc in positions I\ . 

One particularly nice class of positroid cells is the TP or totally positive Schubert cells. 

Definition 2.16. A totally positive Schubert cell is the intersection of a Schubert cell with (Grfc „)>o. 

TP Schubert cells are indexed by J-diagrams such that all boxes are filled with a +. 



The convention of [25] was to place the boundary vertices in clockwise order. 
2 This forbidden pattern is in the shape of a backwards L, and hence is denoted J and pronounced "Le." 

ctually Postnikov's convention was to set ir : (i) = j above, so the decorated permutation we are associating is the 
inverse one to his. 
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3. A Deodhar decomposition of the Grassmannian 

In this section we review Deodhar's decomposition of the flag variety G/B [8], and the parameteriza- 
tions of the components due to Marsh and Rietsch [22] . We then give a Deodhar decomposition of the 
Grassmannian by projecting the usual Deodhar decomposition of the flag variety to the Grassmannian. 

3.1. The flag variety. In this paper we fix the group G = SL„ = SL„(M), a maximal torus T, and 
opposite Borcl subgroups B + and B~ , which consist of the diagonal, upper-triangular, and lower- 
triangular matrices, respectively. We let U + and U~ be the unipotent radicals of B + and B~; these 
are the subgroups of upper-triangular and lower-triangular matrices with l's on the diagonals. For each 
1 < i < ii - 1 we have a homomorphism cf)i : SL 2 — > SL„ such that 



A 



a b 

c d 



\ 



b 
d 



e SL„ 



V 



that is, 4>i replaces a 2 x 2 block of the identity matrix with ^ ^ . Here a is at the (i + l)st 

diagonal entry counting from the southeast corner. 4 We use this to construct 1-parameter subgroups 
in G (landing in U + and U~ , respectively) defined by 



Xi(m) = 4>i 



1 m 
1 



and yi(m) = cf) l 



1 
m 1 



where m G 



Let W denote the Weyl group Nq (T) /T, where Nq (T) is the normalizer of T. The simple reflections 
Si G W are given explicitly by Sj := SiT where s'j := <f>i ( ^ n ) an( i an y w <EW can be expressed 



x 1 

as a product w — s il s i2 . . . s im with m = £(w) factors. We set w = s^s^ . . . Sj m . For G = SL„, we have 
W = & n , the symmetric group on n letters, and Sj is the transposition exchanging i and i + 1. 
We can identify the flag variety G/B with the variety B of Borel subgroups, via 

gB^g-B+ :=gB + g-\ 

We have two opposite Bruhat decompositions of B: 

B= [J B+w-B+= |J B~v-B + . 



wew 



vew 



We define the Richardson variety 



TZ V w := B + w ■ B+ n B~v ■ B+, 



the intersection of opposite Bruhat cells. This intersection is empty unless v < w, in which case it is 
smooth of dimension £(w) — £(v), see [15, 21]. 



Our numbering differs from that in [22] in that the rows of our matrices in SL„ are numbered from the bottom. 
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3.2. Distinguished expressions. We now provide background on distinguished and positive distin- 
guished subexpressions, as in [8] and [22]. We will assume that the reader is familiar with the (strong) 
Bruhat order < on the Weyl group W = & n , and the basics of reduced expressions, as in [4]. 

Let w := s il . . . s im be a reduced expression for w e W. A subexpression v of w is a word obtained 
from the reduced expression w by replacing some of the factors with 1. For example, consider a reduced 
expression in 64, say s 3 S2Sis 3 s 2 s 3 . Then s 3 s 2 1 s 3 s 2 1 is a subexpression of S3S2S1S3S2S3. Given a 
subexpression v, we set v^) to be the product of the leftmost k factors of v, if k > 1, and «( ) = 1- 

Definition 3.1. [22, 8] Given a subexpression v of a reduced expression w = s il s i2 . . . s iml we define 



J' 



= {ke {l,...,m} I U( fe _i) < W( fe )}, 
= {k e {l,...,m} I U( fe _i) = W( fe )}, 
= {fc e {l,...,m} I U( fe _i) > W(fe)}. 



The expression v is called non- decreasing if < uy) for all j = 1, ... , m, e.g. J* = 0. 

Definition 3.2 (Distinguished subexpressions). [8, Definition 2.3] A subexpression v of w is called 
distinguished if we have 

(3.1) V(j) < ^(j-i) s t] for all j e {1, . . . , m}. 

In other words, if right multiplication by Si j decreases the length of vtj-u, then in a distinguished 
subexpression we must have = U(j_i)S^. 

We write v -<! w if v is a distinguished subexpression of w. 

Definition 3.3 (Positive distinguished subexpressions). We call a subexpression v of w a positive dis- 
tinguished subexpression (or a PDS for short) if 

(3.2) < %-i) s »j for a11 3 € {1, . . . ,m}. 
In other words, it is distinguished and non-decreasing. 

Lemma 3.4. [22] Given v < w and a reduced expression w for w, there is a unique PDS v + for v in w. 

3.3. Deodhar components in the flag variety. We now describe the Deodhar decomposition of the 
flag variety. This is a further refinement of the decomposition of G/B into Richardson varieties 1Z V . W . 
Marsh and Rietsch [22] gave explicit paramcterizations for each Deodhar component, identifying each 
one with a subset in the group. 

Definition 3.5. [22, Definition 5.1] Let w = s il . . . s im be a reduced expression for w, and let v be a 
distinguished subexpression. Define a subset G VjW in G by 

U = Xi,{mt)'s~^ if ^ e J', 
mm- ■ ■ Un< ge = Vi e (Pe) if^eJ°, 
ge = s le iileJ°, 

l\ J '\ -> G. 



(3.3) 

There is an obvious map ( 



w = 1 we define G v 



{I}- 



for p{ e I*, mi e R. 

V; w defined by the parameters pi and mi in (3.3). For 



Example 3.6. Let W = 65, w = S2S3S4S1S2S3 and v = s 2 llls 2 l. Then the corresponding element 
g E G VjW is given by g = s 2 y 3 (p2)y4(P3)yi(P4)x2(m 5 )s2 1 y3(P6), which is 



/ 1 











o\ 


P3 


1 














P6 


1 








P2P3 




-m 5 


1 





V 




-P4 





V 
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The following result from [22] gives an explicit paramctrization for the Deodhar component 1Z V ^ 
We will take the description below as the definition of 7£ VjW . 



Proposition 3.7. [22, Proposition 5.2] The map (M*)I J °I x Rl J 'l -> G v , w from Definition 3.5 is an 



7?\ 

isomorphism. The set G v , w lies in U~vr\B + wB + , and the assignment g i-> .g--E> + defines an isomorphism 
(3.4) G v ,w > T^v.-w 

between the subset G v , w of the group, and the Deodhar component 7^ v ,w in the flag variety. 

Suppose that for each w € W we choose a reduced expression w for w. Then it follows from Deodhar's 
work (see [8] and [22, Section 4.4]) that 



(3.5) 



Hv,w = U ^v,w and G/B = |_J f |_J ft v , w J . 

v-<w weW \v^w / 



These are called the Deodhar decompositions oilZ VA1] and G/B. 

Remark 3.8. One may define the Richardson variety 1Z V . W over a finite field F q . In this setting the 
number of points determine the i?-polynomials R v . w (q) = #(7^~u.u)(F 9 )) introduced by Kazhdan and 
Lusztig [14] to give a formula for the Kazhdan-Lusztig polynomials. This was the original motivation for 
Deodhar's work. Therefore the isomorphisms 7£ v ,w = (F*)' Jv x f|/ v ' together with the decomposition 

(3.5) give formulas for the i?-polynomials. 

3.4. Projections of Deodhar components to the Grassmannian. Now we consider the projec- 
tion of the Deodhar decomposition to the Grassmannian Grk, n - Let Wk be the parabolic subgroup 
(si, S2, ■ ■ ■ , Sn-fej ■ • • i s n _i) of W = ©„, and let W k denote the set of minimal- length coset represen- 
tatives of W/Wk- Recall that a descent of a permutation ir is a position j such that ir(j) > ir(j + 1). 
Then W k is the subset of permutations of & n which have at most one descent; and that descent must 
be in position n — k. 

Let nk : G/B — > Grfe >n be the projection from the flag variety to the Grassmannian. For each 
w e W k and v < w, define V VjW = ^k(T^v,w)- Then by work of Lusztig [21], -Kk is an isomorphism on 
V ViW , and we have a decomposition 

(3.6) Gr k , n = U ( U 

wew k \v<w 

Definition 3.9. For each reduced decomposition w for to G W k , and each v -< w, we define the 
(projected) Deodhar component P v ,w = 7rfe(^v,w) C Grk, n - 

Lemma 3.10. [20, Remark 3.12, Lemma 3.13] The decomposition in (3.6) coincides with the positroid 
stratification from Section 2.3. The appropriate bijection between the strata is defined as follows. Let 
Q k denote the set of pairs (v, w) where v <E W, w e W k , and v < w; let Dec^ denote the set of decorated 
permutations in S n with k weak excedances. We consider both sets as partially ordered sets, where 
the cover relation corresponds to containment of closures of the corresponding strata. Then there is 
an order-preserving bijection $ from Q k to Dec* which is defined as follows. Let (u, to) € Q J ■ Then 
<&(v,w) = {it, col) where 7r = t>to -1 . We also let w{v,w) denote $(v,w). To define col, we color any 
fixed point that occurs in one of the positions to(l), w(2), . . . , w(n — k) with the color —1, and color any 
other fixed point with the color 1. 

By Lemma 3.10, "P ViW lies in the positroid stratum SV : • Note that the strata P ViW do not depend on 
the chosen reduced decomposition of w [19, Proposition 4.16]. 
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Now if for each w <E W we choose a reduced decomposition w, then we have 
(3.7) V v , w - |J 7>v, w and Gr k . n = |J ( U Vv <™ ) ' 

v^w weW k \v-!w / 

Proposition 3.7 gives us a concrete way to think about the projected Deodhar components V v 
The projection n k : G/B — > Gr k , n maps each g € G ViW to the span of its leftmost k columns: 

^ 9n,n • • • 9n,n—k+l • • • 9n, ,n—k-\-l • • • Qn.n — k+l\ 



\9i, 



9l,n-k+l 



i y A 



\ 9l,r. 



Alternatively, we may identify A € Gr k ^ n with its image in the Pliicker embedding. Let a denote the 
column vector in R" such that the ith entry from the bottom contains a 1, and all other entries are 0, 
e.g. e n = (1, 0, . . . , 0) T , the transpose of the row vector (1,0,..., 0). Then the projection Tr k maps each 
g G G v , w (identified with g ■ B + e 7£ v , w ) to 

(3.8) 



3 • e„ A . . . A e„_ fe+ i = £ A ju ... tjk (A)e jk A • • • A e n . 

l<ji<-<jk<n 



That is, the Pliicker coordinate Aj 1 ,,,,j k (A) is given by 

A h,-,j k ( A ) = ( e J k A • • • A ejl , 5 • e„ A • • • A e„_ fe+ i), 
where (•, •) is the usual inner product on A fe M™. 

Example 3.11. We continue Example 3.6. Note that w € W k where k = 2. Then the map 7T2 : G v 



G r 2,5 is given by 












( 1 















P3 


1 











9 = 





P6 


1 










P2P3 


Pi - m 5 p 6 


-m 5 


1 







V 


-PaPq 


-P4 





V 



A: 



-P4P6 P2 - ^5P6 P6 1 0\ 
P2P3 p 3 II' 



4. DEODHAR COMPONENTS IN THE GRASSMANNIAN AND GO-DIAGRAMS 

In this section we explain how to index the Deodhar components in the Grassmannian Gr k ^ n by cer- 
tain tableaux called Go-diagrams, which are fillings of Young diagrams by empty boxes, •'s and O's. We 
refer to the symbols • and O as black and white stones. Recall that W k — (si, S2, ■ ■ ■ , Sn-fc, • • • , Sn-i) 
is a parabolic subgroup of W = & n and W k is the set of minimal- length coset representatives of W/W k . 

We fix k and n, and let Q k be the poset whose elements are the boxes in a k x (n — k) rectangle. If 
61 and 62 are two adjacent boxes such that hi is immediately to the left or immediately above b\, we 
have a cover relation b\ < 62 in Q . The partial order on Q k is the transitive closure of <. The middle 
and right diagram in the figure below show two linear extensions (or reading orders) of the poset Q 3 , 
where n = 8. Next we assign a labeling of the boxes of Q k by simple generators Sj, see the left diagram 
of the figure below. If b is a box labelled by s i} we denote the simple generator labeling b by s b :— Sj. 

The following result can be found in [29] . 

Proposition 4.1. Fix k and n. The upper order ideals of Q k are in bijection with elements of W k , and 
the different reading orders of Q k allow us to compute all reduced expressions of elements of W k . More 
specifically let Y be an upper order ideal of Q k , and choose a reading order e for the boxes of Y. Then 
if we read the labels of Y in the order specified by e, we will get a reduced word for some w € W k , and 
this element w does not depend on the choice of e. Therefore we may denote Y by O w . Moreover, if 
we let e vary over all reading orders for boxes of Y, we will obtain all reduced expressions for w. 
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S5 


S 4 


S3 


S2 


Si 


S6 


S5 


S 4 


S3 


S2 


S7 


«6 


S5 


S 4 


S3 



15 


14 


13 


12 


11 


10 


9 


8 


7 


6 


5 


4 


3 


2 


1 



15 


12 


9 


6 


3 


14 


11 


8 


5 


2 


13 


10 


7 


4 


1 



Remark 4.2. The upper order ideals of Q k can be identified with the Young diagrams contained in a 
k x (n — k) rectangle (justified at the upper left), and the reading orders of O w can be identified with 
the reverse standard tableaux of shape O w , i.e. entries decrease from left to right in rows and from top 
to bottom in columns. 

4.1. Go-diagrams and labeled Go-diagrams. Next we will identify distinguished subexpressions of 
reduced words for elements of W k with certain tableaux called Go-diagrams. 

Definition 4.3. Let w be a reduced expression for some w £ W , and consider a distinguished subex- 
pression v of w. Let O w be the Young diagram associated to w and choose the reading order of its 
boxes which corresponds to w. Then for each k G J° we will place a O in the corresponding box of 
O w ; for each k G J* we will place a # in the corresponding box; and for each k € we will leave the 
corresponding box blank. We call the resulting diagram a Go-diagram, and refer to the symbols O and 
# as white and black stones. 

Remark 4.4. Note that a Go-diagram has no black stones if and only if v is a positive distinguished 
subexpression of w. In this case, if we replace the empty boxes by +'s and the white stones by 0's, we 
will get a J-diagram. Therefore, slightly abusing terminology, we will often refer to a Go-diagram with 
no black stones as a J-diagram. 5 Moreover, the Go-diagrams with no black stones are in bijection with 
J-diagrams. See [20, Section 4] for more details. 

Example 4.5. Consider the upper order ideal O w which is Q k itself for 65 and k = 2. Then Q k is the 
poset shown in the left diagram. Let us choose the reading order indicated by the labeling shown in 
the right diagram. 



S3 


«2 


Si 


s 4 


S3 


S2 



6 


5 


4 


3 


2 


1 



Now consider the distinguished subexpressions S2S3S4S1S2S3, IS3S4S1IS3, and IS3IS1IS3 of the re- 
duced expression w = S2S3S4S1S2S3 corresponding to our chosen reading order. Among these three 
subexpressions, the first and second are PDS's. The corresponding Go-diagrams are as follows. 







































• 















The reader might worry that Definition 4.3 has too much dependence on the choice of reduced 
expression w, or cquivalently on the reading order e. However, we have the following result. 

Proposition 4.6. [20, Proposition 4.5] Let D be a Go-diagram in a Young diagram O w . Choose any 
reading order e for the boxes of O w . Let w be the corresponding reduced expression for w. Let v(D) 

5 Since J-diagrams are a special case of Go-diagrams, one might also refer to them as Lego diagrams. 
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be the subexpression of w obtained as follows: if a box b of D contains a black or white stone then 
the corresponding simple generator s b is present in the subexpression, while if b is empty, we omit the 
corresonding simple generator. Then we have the following: 

(1) the element v := v(D) is independent of the choice of reading word e. 

(2) whether v(D) is a PDS depends only on D (and not e). 

(3) whether v(D) is distinguished depends only on D (and not on e). 

Definition 4.7. Let O w be an upper order ideal of Q k , where w G W k and W = S n . Consider a 
Go-diagram D of shape O w ; this is contained in a k x (n — k) rectangle, and the shape O w gives 
rise to a lattice path from the northeast corner to the southwest corner of the rectangle. Label the 
steps of that lattice path from 1 to n; this gives a natural labeling to every row and column of the 
rectangle. We now let v :— v(D), and we define tt'(D) to be the decorated permutation (tt(D),coI) 
where it = ir(D) — vw~ x . The fixed points of ir correspond precisely to rows and columns of the 
rectangle with no +'s. If there are no +'s in the row (respectively, column) labeled by h, then n(h) = h 
and this fixed point gets colored with color 1 (respectively, — 1.) 

Remark 4.8. It follows from Lemma 3.10 that the projected Deodhar component Vd corresponding to 
D is contained in the positroid stratum S n -. ( d) ■ 

Remark 4.9. Recall from Remark 3.8 that the isomorphisms 7\L v ,w — (IFg)'' 7 ° x together with 

the decomposition (3.5) give formulas for the i?-polynomials. Therefore a good characterization of 
Go-diagrams could lead to explicit formulas for the corresponding i?-polynomials. 

If we choose a reading order of O w , then we will also associate to a Go-diagram of shape O w a labeled 
Go-diagram, as defined below. Equivalently, a labeled Go-diagram is associated to a pair (v,w). 

Definition 4.10. Given a reading order of O w and a Go-diagram of shape O w , we obtain a labeled Go- 
diagram by replacing each O with a 1, each • with a —1, and putting a pi in each blank square b, where 
the subscript i corresponds to the label of b inherited from the reading order. 

The labeled Go-diagrams corresponding to the Go-diagrams from Example 4.5 are: 



1 


1 


1 


1 


1 


1 



1 


Pb 


1 


1 


1 


Pi 



-1 


P5 


1 


P3 


1 


Pi 



4.2. Pliicker coordinates for projected Deodhar components. Consider ? v w C Gr^.n, where 
w is a reduced expression for w G W k and v -< w. Here we provide some formulas for the Pliicker 
coordinates of the elements of P ViW , in terms of the parameters used to define G v , w . Some of these 
formulas are related to corresponding formulas for G/B in [22, Section 7]. 

Theorem 4.11. [20, Lemma 5.1, Theorem 5.2] Choose any element A of V VtVr C Grk. n , in other words, 
A = nk(g) for some g G G VjW . Then the lexicographically minimal and maximal nonzero Pliicker 
coordinates of A are Aj and A/< , where 

I = w {n, n — 1, . . . , n — k + 1} and I' = v {n, n — 1, . . . , n — k + 1}. 

Moreover, if we write g = g\ . . . g m as in Definition 3.5, then 

(4.1) A/(A) = (-1) |J *' ]J and A i'( A ) = L 

Note that Aj(A) equals the product of all the labels from the labeled Go-diagram associated to (v, w). 
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This theorem can be extended to provide a formula for some other Plucker coordinates. Let b be 
any box of D. We can choose a linear extension e of the boxes of D which orders the boxes which are 
weakly southeast of b (the inner boxes) before the rest (the outer boxes). This gives rise to a reduced 
expression w and subexpression v, as well as reduced expressions w™ = w ln , v^ 11 = v m , w£ ut = w° 
and v£ ut = v out , which are obtained by restricting w and v to the inner and outer boxes, respectively. 

Theorem 4.12. [20, Theorem 5.6] Let w = Si 1 . . . Si m be a reduced expression for w € W k and v -< w 
and let D be the corresponding Go-diagram. Choose any box b of D, and let v ln = v b n and u> ln = w b n 

Let A = TTk(g) for any g € G v , w , and let / = w{n, n— 1, . . . , n — fc + 1} 



and v out = v° ut and w out 



Define lb — v ln (w ln ) x l E (^')- If we write g — gi ■ ■ ■ g m as in Definition 3.5, then 
(4.2) A 7b (A) = (-l)K-l J] Pi . 



Note that Aj b (A) equals the product of all the labels in the "out" boxes of the labeled Go-diagram. 
Example 4.13. Consider A <G Gr 3 j which is the projection of g e G v , w with 

W = S 3 S 4 S5S 6 S2S3S4S5SlS2S3S4, V = S 3 1 S 5 1 S 2 S 3 1 S 5 1 S 2 S 3 S 4 . 



Then v = s 2 s 4 and n 
follows: 



vw 



(4, 6, 7, 1, 3, 2, 5). The Go-diagram and the labeled Go-diagram are as 



o 


• 


• 




• 




o 


o 




o 




o 



1 


-1 


-1 


P9 


-1 


P7 


1 


1 


Pi 


1 


P2 


1 



467 


456 


245 


234 


167 


156 


125 


123 


127 


125 


125 


123 



Each box b in the diagram on the right show the subset Then we have 

Al, 2 ,3 = -P2P4P7P9, Ai ;2i5 = -piP 7 P9, ^1,2,7 = -P7P9, ^1,5,6 = ~PiP9, 
Ai, 6 ,7 = P9, A 2 , 3 ,4 = -P2P4P7, ^2,4,5 = PiP7 , ^4,5,6 = ~Pi, &lfi,7 = L 

5. SOLITON SOLUTIONS TO THE KP EQUATION AND THEIR CONTOUR PLOTS 

We now explain how to obtain a soliton solution to the KP equation from a point of Grk, n - 

5.1. Prom a point of the Grassmannian to a r- function. We start by fixing real parameters Hi 

p 

such that K\ < k 2 ■ ■ ■ < k„, which are generic, in the sense that the sums ^ K ij are all distinct for 

3 = 1 

any p with 1 < p < n. We also assume that the differences between consecutive Ki's are similar, that 
is, Ki+i — Ki is of order one. 

Let {Ei \ i = 1, . . . , m} be a set of exponential functions in (x, y, t) £ K 3 defined by 

i?i(a;,y,t) := exp 6>i(x, y, where 9i(x,y,t) — k { x + njy + nft. 

Thinking of each Ei as a function of x, we see that the elements {Ei} are linearly independent, because 
their Wronskian determinant with respect to x is non zero: 

Wr(£ 1; ...,E n ):= det[(Sp'- 1) )i< i;j <„] = JJ^- - m) E 1 ■ ■ ■ E n ± 0. 

i<3 

Here E\ i] := d^E t /dx 3 = 4 E i- 

Let A be a full rank k x n matrix. We define a set of functions {/i, . . . , /&} by 

(/i>/2, • ■ • ,fk) T = A ■ (E 1 ,E 2 , ■ ■ ■ ,E n ) T , 
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where (. . ) T denotes the transpose of the vector (...). 

Since the exponential functions {Ei} are linearly independent, we identify them as a basis of MP, and 
then ■ ■ , fk} spans a fc-dimensional subspace. This identification can be seen, more precisely, as 
Ei (1 , Ki, . . . , £ M n . This subspace depends only on which point of the Grassmannian Gr^.n 

the matrix A represents, so we can identify the space of subspaces {/i, . . . , /&} with Grk, n - 

The T-function of A is defined by 

(5-1) r A (x,y,t) :=Wr(/i,/ 2 ,...,/ fc ). 

For I = {i u ...,i k } £ we set 

Ei(x,y,t) :=m(E il ,E i2 ,...,E ik ) = \\{n lm - n it ) E h ■ ■ ■ E ik > 0. 

Km 

n 

Applying the Binet-Cauchy identity to the fact that fj — a jiEi, we get 

(5.2) TA {x,y,t)= £ A I (A)E I (x,y,t). 

It follows that if A e (Gr fe .„)>o, then ta > for all (x,y,t) £ R 3 . 

Thinking of ta as a function of A, we note from (5.2) that the T-function encodes the information 
of the Pliicker embedding. More specifically, if we identify each function Ej with I = . . . ,i k } with 
the wedge product E ix A • • • A E ik (recall the identification Ei o (1, Ki, . . . , k" _1 ) t ), then the map 

t : Grk, n Kp( fc ) _1 , A^ ta has the Pliicker coordinates as coefficients. 

5.2. From the T-function to solutions of the KP equation. The KP equation for u(x,y,t) 

d / ^du du d 3 u\ 2,® 2u 
dx \ dt dx dx 3 J dy 2 

was proposed by Kadomtsev and Petviashvili in 1970 [13], in order to study the stability of the soliton 
solutions of the Korteweg-de Vries (KdV) equation under the influence of weak transverse perturbations. 
The KP equation can be also used to describe two-dimensional shallow water wave phenomena (sec for 
example [17]). This equation is now considered to be a prototype of an intcgrablc nonlinear partial 
differential equation. For more background, see [24, 9, 1, 12, 23]. 

It is well known (see [12, 5, 6, 7]) that the T-function defined in (5.1) provides a soliton solution of 
the KP equation, 

d 2 

(5.3) u A (x,y,t) = 2^ ln t a (x, y,t). 

It is easy to show that if A £ (Gr^„)> , then such a solution ua(x, y, t) is regular for all (a;, y, t) £ R 3 . 
For this reason we are interested in solutions UA(x,y,t) of the KP equation which come from points A 
of (Grk, n )>o- Throughout this paper when we speak of a soliton solution to the KP equation, we will 
mean a solution UA{x,y,t) which has form (5.3). 

5.3. Contour plots of soliton solutions. One can visualize a solution ua{x, y, t) to the KP equation 
by drawing level sets of the solution in the a;y-plane, when the coordinate t is fixed. For each r£l, 
we denote the corresponding level set by 

C r (t) := {(x,y) £ W 2 : u A (x,y,t) = r}. 

Figure 2 depicts both a three-dimensional image of a solution Ua{x, y,t), as well as multiple level sets 
C r . These level sets are lines parallel to the line of the wave peak. 
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Figure 2. A line-soliton solution ua{x,u, t) where A — (1,1) € (Gri,2)>0i depicted 
via the 3-dimensional profile UA(x,y,t), and the level sets of UA(x,y,t) for some t. 
represents the dominant exponential in each region. 

To study the behavior of UA(x,y,t) for A £ Sm C Gr^ „, we consider the dominant exponentials in 
the r-function (5.2) at each point {x,y,t). First we write the r-function in the form 

T A (x,y,t)= £ Aj(A)Ej(x,y,t) 

= E exp ( T,(K ji x + K 2 j .y + K 3 j .t) +ln(Aj(A)Kj) 

JeM \i=l 

where Kj :— Yie<m( K j m ~ K je) > 0- Note that in general the terms ln(Aj(A)Kj) could be imaginary 

when some Aj(A) are negative. 

Since we are interested in the behavior of the soliton solutions when the variables (x, y, t) are on a 

large scale, we rescale the variables with a small positive number e, 

x y t 

x — ► -, y — > -, t — > -. 
e e e 

This leads to 

T%(x,y,t)= expf- 'E(Kj i x + K 2 ji y + nit) + ln{A j{A)Kj) 

JEM \ e i=l 

Then we define a function /a{x, y, t) as the limit 

f A (x,y,t) = Umeln(r|(x, y,i)) 

(5.4) r fe i 

= max < V (/s,-, a; + k^.w + > . 
JeM J Jz 1% J 

Since the above function depends only on the collection M., we also denote it as fM(x,y,t). 

Definition 5.1. Given a solution UA{x,y,t) of the KP equation as in (5.3), we define its contour plot 
C{ua) to be the locus in M 3 where /a(x, y, t) is not linear. If we fix t — t 0l then we let C to (uA) be the 
locus in K 2 where Ja(x, y,t — to) is not linear, and we also refer to this as a contour plot. Because these 
contour plots depend only on M. and not on A, we also refer to them as C{M) and C to (M). 

Remark 5.2. The contour plot approximates the locus where \uA(x,y,t)\ takes on its maximum values 
or is singular. 

Remark 5.3. Note that the contour plot generated by the function /a (a;, y,t) at t — consists of a set 
of semi-infinite lines attached to the origin (0, 0) in the xy-plane. And if t\ and £2 have the same sign, 
then the corresponding contour plots C tl (M) and C t2 (M) are self-similar. 
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Also note that because our definition of the contour plot ignores the constant terms \n(Aj(A)Kj), 
there are no phase-shifts in our picture, and the contour plot for /a{x, y, t) = ,/>((#, y, t) does not 
depend on the signs of the Pliicker coordinates. 

It follows from Definition 5.1 that C{ua) and C ta {uA) are piecewise linear subsets of R 3 and R 2 , 
respectively, of codimension 1. In fact it is easy to verify the following. 

Proposition 5.4. [19, Proposition 4.3] If each Ki is an integer, then C(ua) is a tropical hypersurface in 
R 3 , and C ta {uA) is a tropical hypersurface (i.e. a tropical curve) in R 2 . 

The contour plot Ct {uA) consists of line segments called line-solitons, some of which have finite 
length, while others are unbounded and extend in the y direction to ±oo. Each region of the complement 
of Ct (ua) in R 2 is a domain of linearity for Ja(x, y, t), and hence each region is naturally associated to 
a dominant exponential Aj(A)Ej(x,y,t) from the r-function (5.2). We label this region by J or Ej. 
Each line-soliton represents a balance between two dominant exponentials in the r-function. 

Because of the genericity of the K-parameters, the following lemma is immediate. 

Lemma 5.5. [7, Proposition 5] The index sets of the dominant exponentials of the r-function in adjacent 
regions of the contour plot in the xy-plane are of the form {i, l 2 , ■ ■ ■ , /&} and {j, l 2 , ■ ■ ■ , h}- 

We call the line-soliton separating the two dominant exponentials in Lemma 5.5 a line-soliton of type 
Its equation is 

(5.5) X + (Ki + Kj)y + (k 2 + KiKj + K 2 )t = 0. 

Remark 5.6. Consider a line-soliton given by (5.5). Compute the angle ^[ij] between the positive y-axis 
and the line-soliton of type measured in the counterclockwise direction, so that the negative cc-axis 
has an angle of ^ and the positive x-axis has an angle of — |. Then tan^^ j] = Ki + Kj, so we refer to 
the slope of the line-soliton (see Figure 2). 

In Section 7 we will explore the combinatorial structure of contour plots, that is, the ways in which 
line-solitons may interact. Generically we expect a point at which several line-solitons meet to have 
degree 3; we regard such a point as a trivalent vertex. Three line-solitons meeting at a trivalent vertex 
exhibit a resonant interaction (this corresponds to the balancing condition for a tropical curve). See 
[19, Section 4.2]. One may also have two line-solitons which cross over each other, forming an X-shape: 
we call this an X-crossing, but do not regard it as a vertex. See Figure 4. Vertices of degree greater 
than 4 are also possible. 

Definition 5.7. Let i < j < k < £ be positive integers. An A-crossing involving two line-solitons of 
types [i, k] and [j, l\ is called a black X-crossing. An A-crossing involving two line-solitons of types [i, j] 
and [k,£], or of types [i,£] and [j, k], is called a white X-crossing. 

Definition 5.8. A contour plot C((ua) is called generic if all interactions of line-solitons are at trivalent 
vertices or are A-crossings. 

6. Unbounded line-solitons at y > and i/«0 

In this section we explain that the unbounded line-solitons at \y\ 3> of a contour plot Ct(uA) are 
determined by which positroid stratum contains A. Conversely, the unbounded line-solitons of Ct(uA) 
determine which positroid stratum A lies in. 

Theorem 6.1. [20, Theorem 8.1] Let A e Grk, n he in the positroid stratum S^., where ir' — (n,col). 
Consider the contour plot Ct{uA) for any time t. Then the excedances (respectively, nonexcedances) 
of 7r are in bijection with the unbounded line-solitons of C t {uA) at y ^> (respectively, y <^ 0). More 
specifically, in C t (uA), 
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(a) there is an unbounded line-soliton of [i, h]-type at y ^ if and only if w(i) = h for i < h, 

(b) there is an unbounded line-soliton of [i, h]-type at y <C if and only if ir(h) — i for i < h. 

Therefore w determines the unbounded line-solitons at y 3> and y <C of C t (uA) for any time t. 

Conversely, given a contour plot Ct{uA) at any time t where A G Grk, n , one can construct ir' = (ir, col) 
such that A G S n - as follows. The excedances and nonexcedances of tt are constructed as above from 
the unbounded line-solitons. If there is an h G [n] such that h G J for every dominant exponential Ej 
labeling the contour plot, then set n(h) = h with col(h) = 1. If there is an h G [n] such that h £ J for 
any dominant exponential Ej labeling the contour plot, then set ir(h) — h with col(h) = —1. 

Remark 6.2. Chakravarty and Kodama [5, Prop. 2.6 and 2.9] and [7, Theorem 5] associated a derange- 
ment to each irreducible element A in the totally non-negative part (Grk,n)>o of the Grassmannian. 
Theorem 6.1 generalizes their result by dropping the hypothesis of irreducibility and extending the 
setting from (Gr fe! „)> to Gr k , n . 

Example 6.3. Consider some A G Gr 4i9 which is the projection of an element g G G v , w with 

W = S- ! S^S4 L S^ ) SQS'jS2SiS^S(,SiS2S^S4 L Sc ) and V = S 7 llS5lS7S2lS4lllS2lS4S5- 

Then v = 1 and n = vw^ 1 = (6, 7, 1, 8, 2, 3, 9, 4, 5). The matrix g G G v , w is given by 

9 = s 7 y8(P2)yi(P3)s5y6(P5)x 7 (m 6 )s^ 1 S2y 3 (Ps)s4y5(Pw)y6(Pii) 
■ yi(pi2)x2(m 13 )s2 1 y 3 (p 1 4)x 4 (m 15 )s4 1 x 5 (m 16 )s5 1 . 
The Go-diagram and the labeled Go-diagram are as follows: 



• 


• 




• 








o 




o 


• 




o 








o 









-1 


-1 


Pl4 


-1 


Pl2 


Pll 


PlO 


1 


P8 


1 


-1 


P5 


1 


P3 




P2 


1 









The A-matrix is then given by 

/-P12P14 







913 

PsPioPu 





Pl4 







915 

Pll(.P3 +Pw) 
-P3P5 





-77116 
Pll 
~P5 




1 

1 

— 77l 6 

P2 



\ 




1/ 



where the matrix entry q 13 = -mi 3 pu + mi 5 p s - m 16 p 8 pio and 1715 = m 15 - m le (p 3 +p w ). In Figure 3, 
we show contour plots C((ua) for the solution ua(x, y, t) at t = —10, 0, 10, using the choice of parameters 
(«i, . . . , Kg) = (—5, —3, —2, —1, 0, 1, 2, 3, 4), pj = 1 for all j, and to/ = for all I. Note that: 

(a) For y ^> 0, there are four unbounded line-solitons, whose types from right to left are: 

[1,6], [2,7], [4,8], and [7,9] 

(b) For y <C 0, there are five unbounded line-solitons, whose types from left to right are: 

[1,3], [2,5], [3,6], [4,8], and [5,9] 

We can see from this example that the line-solitons for y >• correspond to the excedances in n = 
(6, 7, 1, 8, 2, 3, 9, 4, 5), while those for i/«0 correspond to the nonexcedances. 

Note that if there are two adjacent regions of the contour plot whose Plucker coordinates have different 
signs, then the line-soliton separating them is singular. For example, the line-soliton of type [4, 8] (the 
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t = - 10 t = t= 10 



Figure 3. Example of contour plots Ct{uA) for A G GV49. The contour plots are 
obtained by "Plot3D" of Mathematica (see the details in the text). 

second soliton from the left in y ^ 0) is singular, because the Pliicker coordinates corresponding to the 
(dominant exponentials of the) adjacent regions are 

Ai, 2 ,4,9 = P3P5PSP10P11P12PH = 1 and Ai, 2 ,8,9 = -P8P10P11P12PU = -1. 
7. Soliton graphs and generalized plabic graphs 

The following notion of soliton graph forgets the metric data of the contour plot, but preserves the 
data of how line-solitons interact and which exponentials dominate. 

Definition 7.1. Let A G Gr^ n and consider a generic contour plot C t {uA) for some time t. Color a 
trivalent vertex black (respectively, white) if it has a unique edge extending downwards (respectively, 
upwards) from it. We preserve the labeling of regions and edges that was used in the contour plot: we 
label a region by Ei if the dominant exponential in that region is AiEj, and label each line-soliton 
by its type (see Lemma 5.5). We also preserve the topology of the graph, but forget the metric 
structure. We call this labeled graph with bicolored vertices the soliton graph G to (uA)- 

Example 7.2. We continue Example 6.3. Figure 4 contains the same contour plot Ct(uA) as that at the 
left of Figure 3. One may use Lemma 5.5 to label all regions and edges in the soliton graph. After 
computing the Pliicker coordinates, one can identify the singular solitons, which are indicated by the 
dotted lines in the soliton graph. 

We now describe how to pass from a soliton graph to a generalized plabic graph. 

Definition 7.3. A generalized plabic graph is an undirected graph G drawn inside a disk with n boundary 
vertices labeled {1, . . . , n}. We require that each boundary vertex i is either isolated (in which case it is 
colored with color 1 or —1), or is incident to a single edge; and each internal vertex is colored black or 
white. Edges are allowed to cross each other in an X-crossing (which is not considered to be a vertex). 

By Theorem 6.1, the following construction is well-defined. 

Definition 7.4. Fix a positroid stratum S^- of Grk^n where w — (tt,coI). To each soliton graph C 
coming from a point of that stratum we associate a generalized plabic graph Pl{C) by: 

• embedding C into a disk, so that each unbounded line-soliton of C ends at a boundary vertex] 

• labeling the boundary vertex incident to the edge with labels i and 7r(i) by 7r(i); 

• adding an isolated boundary vertex labeled h with color 1 (respectively, —1) whenever h G J 
for each region label Ej (respectively, whenever h ^ J for any region label Ej); 
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Figure 4. Example of a contour plot Ct(uA), its soliton graph C = Gt(v,A), arid its 
generalized plabic graph Pl{C). The parameters used are those from Example 6.3. In 
particular, . . . , Kg) = (—5, —3, —2, —1, 0, 1, 2, 3, 4), and 7r = (6, 7, 1, 8, 2, 3, 9, 4, 5). 

• forgetting the labels of all edges and regions. 

See Figure 4 for a soliton graph C together with the corresponding generalized plabic graph Pl{C). 

Definition 7.5. Given a generalized plabic graph G, the trip T is the directed path which starts at the 
boundary vertex i, and follows the "rules of the road": it turns right at a black vertex, left at a white 
vertex, and goes straight through the X-crossings. Note that Tj will also end at a boundary vertex. If 
i is an isolated vertex, then T starts and ends at i. Define 7Tg(«) = j whenever T ends at j. It is not 
hard to show that ttq is a permutation, which we call the trip permutation. 

We use the trips to label the edges and regions of each generalized plabic graph. 

Definition 7.6. Given a generalized plabic graph G, start at each non-isolated boundary vertex i and 
label every edge along trip Tj with i. Such a trip divides the disk containing G into two parts: the part 
to the left of T, and the part to the right. Place an i in every region which is to the left of Tj. If h is 
an isolated boundary vertex with color 1, put an h in every region of G. After repeating this procedure 
for each boundary vertex, each edge will be labeled by up to two numbers (between 1 and n) , and each 
region will be labeled by a collection of numbers. Two regions separated by an edge labeled by both i 
and j will have region labels S and (S \ {i}) U {j}. When an edge is labeled by two numbers i < j, we 
write [i,j] on that edge, or {i, j} or {j, i} if we do not wish to specify the order of i and j. 

Theorem 7.7. [19, Theorem 7.6] Consider a soliton graph C = G t {uA) coming from a point A of a 
positroid stratum S K --, where 7T — (n,col). Then the trip permutation of Pl{C) is tt, and by labeling 
edges of Pl(C) according to Definition 7.6, we will recover the original edge and region labels in C. 

We invite the reader to verify Theorem 7.7 for the graphs in Figure 4. 

Remark 7.8. By Theorem 7.7, we can identify each soliton graph C with its generalized plabic graph 
Pl(C). From now on, we will often ignore the labels of edges and regions of a soliton graph, and simply 
record the labels on boundary vertices. 

8. The contour plot for t <c 

Consider a matroid stratum Sm contained in the Deodhar component Sd, where D is the corre- 
sponding or Go-diagram. From Definition 5.1 it is clear that the contour plot associated to any A E Sm 
depends only on Ai, not on A. In fact for t < a stronger statement is true - the contour plot for 
any A <E Sm C Sd depends only on D, and not on M.. In this section we will explain how to use D to 
construct first a generalized plabic graph G_(D), and then the contour plot C t (M) for t ^ 0. 
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8.1. Definition of the contour plot for t <C 0. Recall from (5.4) the definition of fM(x,y,t). To 
understand how it behaves for t <C 0, let us rescale everything by t. Define x = | and y = | , and set 



that is, KiX + n 2 y + n\t = t(f>i(x, y). Note that because t is negative, x and y have the opposite signs of 
x and y. This leads to the following definition of the contour plot for t <C 0. 

Definition 8.1. We define the contour plot C-oo(M) to be the locus in R 2 where 



Remark 8.2. After a 180° rotation, C- ao (A4) is the limit of C t (uA) as t — > — oo, for any ^4 e S^f. 
Note that the rotation is required because the positive x-axis (respectively, y-axis) corresponds to the 
negative x-axis (respectively, y-axis). 

Definition 8.3. Define Vi.e, m to be the point in R 2 where <pi(x,y) = (j)g(x,y) = <p m (x,y). A simple 
calculation yields that the point Vi t t, m nas the following coordinates in the xy-plane: 



Some of the points Wj^ lTO € R 2 correspond to trivalent vertices in the contour plots we construct; 
such a point is the location of the resonant interaction of three line-solitons of types [£,m] and 

[i, m] (see Theorem 8.5 below). Because of our assumption on the genericity of the K-parameters, those 
points are all distinct. 

8.2. Main results on the contour plot for (<0. This section is based on [20, Section 10], and 
generalizes the results of [19, Section 8] to a soliton solution coming from an arbitrary point of the 
real Grassmannian (not just the non- negative part). We start by giving an algorithm to construct a 
generalized plabic graph G-(D), which will be used to construct C- oc (A4). Figure 5 illustrates the 
steps of Algorithm 8.4, starting from the Go-diagram of the Deodhar component Sd where D is as in 
the upper left corner of Figure 5. 

Algorithm 8.4. [20, Algorithm 10.4] From a Go-diagram D to G_{D): 

(1) Start with a Go-diagram D contained in a k x (n — k) rectangle, and replace each O, •, and 
blank box by a cross, a cross, and a pair of elbows, respectively. Label the n edges along the 
southeast border of the Young diagram by the numbers 1 to n, from northeast to southwest. 
The configuration of crosses and elbows forms n "pipes" which travel from the southeast border 
to the northwest border; label the endpoint of each pipe by the label of its starting point. 

(2) Add a pair of black and white vertices to each pair of elbows, and connect them by an edge, as 
shown in the upper right of Figure 5. Forget the labels of the southeast border. If there is an 
endpoint of a pipe on the east or south border whose pipe starts by going straight, then erase 
the straight portion preceding the first elbow. If there is a horizontal (respectively, vertical) 
pipe starting at i with no elbows, then erase it, and add an isolated boundary vertex labeled i 
with color 1 (respectively, —1). 

(3) Forget any degree 2 vertices, and forget any edges of the graph which end at the southeast 
border of the diagram. Denote the resulting graph G-(D). 

(4) After embedding the graph in a disk with n boundary vertices (including isolated vertices) we 
obtain a generalized plabic graph, which we also denote G_(D). If desired, stretch and rotate 
G-(D) so that the boundary vertices at the west side of the diagram are at the north instead. 



4>i{x,y) = k z x + K 2 y + k*, 




is not linear . 



(tviKg + KiK m + KgK m7 ^Ki -(- Kg ~\- ft m )). 



The following is the main result of this section. 
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Figure 5. Construction of the generalized plabic graph G-(D) associated to the Go- 
diagram D. The labels of the regions of the graph indicate the index sets of the 
corresponding Pliicker coordinates. Using the notation of Definition 4.7, we have 
tt(D) = vw- 1 = (5, 7, 1, 6, 8, 3, 4, 2). 



Theorem 8.5. [20, Theorem 10.6] Choose a matroid stratum Sm an d let Sd be the Deodhar component 
containing Sm- Recall the definition of tt(D) from Definition 4.7. Use Algorithm 8.4 to obtain G-(D). 
Then G-(D) has trip permutation tt(D), and we can use it to explicitly construct C-oo(M) as follows. 
Label the edges of G_ (D) according to the rules of the road. Label by Vi£ tTn each trivalent vertex 
which is incident to edges labeled [i,£\, [i,rn], and [i,m], and give that vertex the coordinates (x,y) = 
(KiKe + KiK m + K£K m , — («i + K( + K m ))- Replace each edge labeled which ends at a boundary 
vertex by an unbounded line-soliton with slope Ki + Kj. (Each edge labeled between two trivalent 
vertices will automatically have slope «i + Kj.) In particular, C~ 00 {M) is determined by D. Recall from 
Remark 8.2 that after a 180° rotation, C_ 00 (.M) is the limit of Ct{uA) as t — > — oo, for any A £ Sm- 

Remark 8.6. Since the contour plot C_ 00 (A / J) depends only on D, we also refer to it as C_oo(-D)- 

Remark 8.7. The results of this section may be extended to the case t 3> by duality considerations 
(similar to the way in which our previous paper [19] described contour plots for both t and t>0). 
Note that the Deodhar decomposition of Grk, n depends on a choice of ordered basis (ei, . . . , e„). Using 
the ordered basis (e„, . . . , e\) instead and the corresponding Deodhar decomposition, one may explicitly 
describe contour plots at t 3> 0. 

Remark 8.8. Depending on the choice of the parameters Ki, the contour plot C_oo(D) may have a 
slightly different topological structure than the soliton graph G_(Z?). While the incidences of line- 
solitons with trivalent vertices are determined by G-(D), the locations of X-crossings may vary based 
on the Ki's. More specifically, changing the «j's may change the contour plot via a sequence of slides, 
see [20, Section 11]. 
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8.3. A-crossings in the contour plots. Recall the notions of black and white A-crossings from 
Definition 5.7. In [19, Theorem 9.1], we proved that the presence of A-crossings in contour plots at 
\t\ » implies that there is a two-term Pliicker relation. 

Theorem 8.9. [19, Theorem 9.1] Suppose that there is an A-crossing in a contour plot Ct(v,A) for some 
A e Grk.n where \t\ ^S> 0. Let I\, I2, h, and I4 be the /c-element subsets of {1, ... , n} corresponding to 
the dominant exponentials incident to the A-crossing listed in circular order. 

• If the A-crossing is white, we have A Il (A)A h (A) = A l2 {A)A Ii {A). 

• If the A-crossing is black, we have A I± (A) A/ 3 (A) = -A l2 (A)A Ii (A). 

The following corollary is immediate. 

Corollary 8.10. If there is a black A-crossing in a contour plot at t < or t > 0, then among the 
Pliicker coordinates associated to the dominant exponentials incident to that black A-crossing, three 
must be positive and one negative, or vice-versa. 

Corollary 8.11. Let D be a J-diagram, that is, a Go-diagram with no black stones. Let A £ Sd and 
t <C 0. Choose any K\ < ■ ■ ■ < n n . Then the contour plot C t {uA) can have only white A-crossings. 

9. Total positivity, regularity, and cluster algebras 

In this paper we have been studying solutions ua(x, y, t) to the KP equation coming from points of the 
real Grassmannian. Among these solutions, those coming from the totally non-negative part of Grk.n 
are especially nice. In particular, (Grk, n )>o parameterizes the regular soliton solutions coming from 
Gr/;. „, see Theorem 9.1. This result provides an important motivation for studying the soliton solutions 
coming from (Grk, n )>o- After discussing the regularity result below, we will discuss positivity tests for 
elements of the Grassmannian, as well as the connection between soliton solutions from (Grk, n )>o and 
the theory of cluster algebras. 

Theorem 9.1. [20, Theorem 12.1] Fix parameters K\ < ■ ■ ■ < K n and an element A e Grk.n- Consider 
the corresponding soliton solution UA{x,y,t) of the KP equation. This solution is regular at t <C if 
and only if A € (Grk, n )>o- Therefore this solution is regular for all times t if and only if A € (Grfc ; „)>o- 

If one is studying total positivity on the Grassmannian, then a natural question is the following: given 
A e Gr fei „, how many Pliicker coordinates, and which ones, must one compute, in order to determine 
that A e (Gr^„)> ? This leads to the following notion of positivity test. 

Definition 9.2. Consider the Deodhar component Sjj C Grj. „, where I? is a Go-diagram. A collection 
J of fc-element subsets of {1, 2, . . . , n} is called a positivity test for Sd if for any A e So, the condition 
that A/(A) > for all / e J implies that A e {Gr k . n )> . 

It turns out that the collection of Pliicker coordinates corresponding to dominant exponentials in 
contour plots at t <C provide positivity tests for positroid strata. 

Theorem 9.3. [20, Theorem 12.9] Let A £ Sd C Grk, n , where D is a J-diagram, and let t « 0. If 
A j (A) > for each dominant exponential Ej in the contour plot Ct(v,A), then A € (Grk. n )>o- In 
other words, the Pliicker coordinates corresponding to the dominant exponentials in C t {uA) comprise a 
positivity test for Sd- 

Remark 9.4. Recall that the Deodhar components Sd have non-empty intersection with (Grk, n )>o 
unless D is a J-diagram. Therefore Theorem 9.3 restricts to the case that D is a J-diagram. 



22 



YUJI KODAMA AND LAUREN WILLIAMS 



9.1. TP Schubert cells and reduced plabic graphs. In this section we provide a new characteri- 
zation of the so-called reduced plabic graphs [25, Section 12]. This will allow us to make a connection 
to cluster algebras in Section 9.2. 

We will always assume that a plabic graph is leafless, i.e. that it has no non-boundary leaves, and 
that it has no isolated components. In order to define reduced, we first define some local transformations 
of plabic graphs. 

(Ml) SQUARE MOVE. If a plabic graph has a square formed by four trivalent vertices whose colors 
alternate, then we can switch the colors of these four vertices. 

(M2) UNICOLORED EDGE CONTRACTION /UNCONTRACTION. If a plabic graph contains an 
edge with two vertices of the same color, then we can contract this edge into a single vertex with the 
same color. We can also uncontract a vertex into an edge with vertices of the same color. 

(M3) MIDDLE VERTEX INSERTION/REMOVAL. If a plabic graph contains a vertex of degree 2, 
then we can remove this vertex and glue the incident edges together; on the other hand, we can always 
insert a vertex (of any color) in the middle of any edge. 




FIGURE 6. A square move; a unicolored edge contraction; a middle vertex insertion/ removal 

(Rl) PARALLEL EDGE REDUCTION. If a network contains two trivalent vertices of different 
colors connected by a pair of parallel edges, then we can remove these vertices and edges, and glue the 
remaining pair of edges together. 



Figure 7. Parallel edge reduction 

Definition 9.5. [25] Two plabic graphs are called move-equivalent if they can be obtained from each 
other by moves (M1)-(M3). The move-equivalence class of a given plabic graph G is the set of all plabic 
graphs which are move-equivalent to G. A leafless plabic graph without isolated components is called 
reduced if there is no graph in its move-equivalence class to which we can apply (Rl). 

Theorem 9.6. [25, Theorem 13.4] Two reduced plabic graphs which each have n boundary vertices are 
move-equivalent if and only if they have the same trip permutation. 

Our new characterization of reduced plabic graphs is as follows. 

Definition 9.7. We say that a (generalized) plabic graph has the resonance property, if after labeling 
edges via Definition 7.6, the set E of edges incident to a given vertex has the following property: 

• there exist numbers ii < i 2 < ■ ■ ■ < i m such that when we read the labels of E, we see the 
labels [ii, i 2 ], [« 2 , 13], • ■ • , [i m _i,i m ], [«i,i TO ] appear in counterclockwise order. 

Theorem 9.8. [19, Theorem 10.5]. A plabic graph is reduced if and only if it has the resonance property. 6 

Corollary 9.9. [19, Corollary 10.9] Suppose that A lies in a TP Schubert cell, and that for some time 
t, G t {uA) is generic with no X-crossings. Then G t {uA) is a reduced plabic graph. 

^Recall from Definition 2.11 that our convention is to label boundary vertices of a plabic graph 1, 2, . . . , n in counter- 
clockwise order. If one chooses the opposite convention, then one must replace the word counterclockwise in Definition 
9.7 by clockwise. 
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9.2. The connection to cluster algebras. Cluster algebras are a class of commutative rings, intro- 
duced by Fomin and Zelevinsky [10], which have a remarkable combinatorial structure. Many coordinate 
rings of homogeneous spaces have a cluster algebra structure: as shown by Scott [28], the Grassmannian 
is one such example. 

Theorem 9.10. [28] The coordinate ring of (the affine cone over) Gr^^ n has a cluster algebra structure. 
Moreover, the set of Plucker coordinates whose indices come from the labels of the regions of a reduced 
plabic graph for (Grk, n )>o comprises a cluster for this cluster algebra. 

Remark 9.11. In fact [28] used the combinatorics of alternating strand diagrams, not reduced plabic 
graphs, to describe clusters. However, alternating strand diagrams are in bijection with reduced plabic 
graphs [25]. 

Theorem 9.12. [19, 10.12] The set of Plucker coordinates labeling regions of a generic trivalent soliton 
graph for the TP Grassmannian is a cluster for the cluster algebra associated to the Grassmannian. 

Conjecturally, every positroid cell <5>£"™ of the totally non-negative Grassmannian also carries a 
cluster algebra structure, and the Plucker coordinates labeling the regions of any reduced plabic graph 
for S^ nn should be a cluster for that cluster algebra. In particular, the TP Schubert cells should 
carry cluster algebra structures. Therefore we conjecture that Theorem 9.12 holds with "Schubert 
cell" replacing "Grassmannian." Finally, there should be a suitable generalization of Theorem 9.12 for 
arbitrary positroid cells. 

9.3. Soliton graphs for (Gr2,„)>o and triangulations of an rt-gon. In this section we must use a 
slightly more general setup for r-functions TA{t\, . . . , t m ) and soliton solutions UA{t\, . . . , t m ), in which 
the exponential functions Ei and Ej are functions of variables t\,. . . ,t m , and m > n. Here t\ = x, 
h = V, h = t, and the other t^s are referred to as higher times. See [19, Section 3] for details. 



1 6 




3 4 



Figure 8. Algorithm 9.13, starting from a triangulation of a hexagon. The right figure 
shows the corresponding contour plot for the soliton solution for (Gr2,e)>o- 

Algorithm 9.13. [19, Algorithm 12.1] Let T be a triangulation of an n-gon P, whose n vertices are labeled 
by the numbers 1, 2, . . . , n, in counterclockwise order. Therefore each edge of P and each diagonal of 
T is specified by a pair of distinct integers between 1 and n. The following procedure yields a labeled 
graph *(T). 

(1) Put a black vertex in the interior of each triangle in T. 
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(2) Put a white vertex at each of the n vertices of P which is incident to a diagonal of T; put a 
black vertex at the remaining vertices of P. 

(3) Connect each vertex which is inside a triangle of T to the three vertices of that triangle. 

(4) Erase the edges of T, and contract every pair of adjacent vertices which have the same color. 
This produces a new graph G with n boundary vertices, in bijection with the vertices of the 
original n-gon P. 

(5) Add one unbounded ray to each of the boundary vertices of G, so as to produce a new (planar) 
graph ty(T). Note that ^f(T) divides the plane into regions; the bounded regions correspond to 
the diagonals of T, and the unbounded regions correspond to the edges of P. 




See Figure 8. Our main result in this section is the following. 

Theorem 9.14. [19, Theorem 12.1] The graphs W(T) constructed above are soliton graphs for (Gr2, n )>o, 
and conversely, up to (M2)-equivalence, any trivalent generic soliton graph for (Gr2, n )>o comes from 
this construction. Moreover, one can realize each graph $>(T) by either: 

• choosing an arbitrary A £ (Gr 2i n)>o an d varying the higher times t^, .. . ,t n appropriately, or 

• fixing an arbitrary collection of higher times ta, . . . , t n , and using the torus action to choose an 
appropriate A £ (Gr2, n )>o- 

Figure 9 shows the five triangulations of a pentagon, together with (contour plots corresponding to) 
the five different soliton graphs which one may obtain from (Gr2^)>o- 

Remark 9.15. The process of flipping a diagonal in the triangulation corresponds to a mutation in the 
cluster algebra. In the terminology of reduced plabic graphs, a mutation corresponds to the square 
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move (Ml). In the setting of KP solitons, each mutation may be considered as an evolution along a 
particular flow of the KP hierarchy defined by the symmetries of the KP equation. 

Remark 9.16. It is known already that the set of reduced plabic graphs for the TP part of Gr 2 . n all 
have the form given by Algorithm 9.13. And by Corollary 9.9, every generic soliton graph is a reduced 
plabic graph. Therefore it follows immediately that every soliton graph for the TP part of Gr 2n must 
have the form of Algorithm 9.13. To prove Theorem 9.14, one must also show that every outcome of 
Algorithm 9.13 can be realized as a soliton graph. 



10. THE INVERSE PROBLEM FOR SOLITON GRAPHS 

The inverse problem for soliton solutions of the KP equation is the following: given a time t together 
with the contour plot Ct(uA) of a soliton solution, can one reconstruct the point A of Gru y n which gave 
rise to the solution? Note that solving for A is desirable, because this information would allow us to 
compute the entire past and the entire future of the soliton solution. 

Using the cluster algebra structure for Grassmannians, we have the following. 

Theorem 10.1. [19, Theorem 11.2] Consider a generic contour plot Ct(uA) of a soliton solution which has 
no A-crossings, and which comes from a point A of the totally positive Grassmannian at an arbitrary 
time t. Then from the contour plot together with t we can uniquely reconstruct the point A. 

Using the description of contour plots of soliton solutions when t <C 0, we have the following. 

Theorem 10.2. [19, Theorem 11.3] Fix K\ < ■ ■ ■ < K n as usual. Consider a generic contour plot C t {uA) 
of a soliton solution coming from a point A of a positroid cell S^ n , for t <C 0. Then from the contour 
plot together with t we can uniquely reconstruct the point A. 

10.1. Non-uniqueness of the evolution of the contour plots for t >• 0. In contrast to the totally 
non-negative case, where the soliton solution can be uniquely determined by the information in the 
contour plot at t -C 0, if we consider arbitrary points A G Gr^_ ni we cannot solve the inverse problem. 

Consider ieSflC Gr^, n - If the contour plot C-oo(D) is topologically identical to G-(D), then the 
contour plot has almost no dependence on the parameters mj from the parameterization of Sd- This 
is because the Pliicker coordinates corresponding to the regions of C_oo(D) (representing the dominant 
exponentials) are either monomials in theft's (see [19, Section 5] and [19, Remark 10.5]), or determined 
from these by a "two-term" Pliicker relation. 

Therefore it is possible to choose two different points A and A' in Sd C Grk, n whose contour plots 
for a fixed m < . . . n n and fixed t <C are identical (up to some exponentially small difference); we use 
the same parameters Pi but different parameters rrij for defining A and A' . However, as t increases, 
those contour plots may evolve to give different patterns. 

Consider the Deodhar stratum Sd C Gr 2i , corresponding to 

w = s 2 s 3 sis 2 and v = s 2 lls 2 . 
The Go-diagram and labeled Go-diagram are given by 



• 
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The matrix g is calculated as g — s 2 y^{p 2 )yi{pz)x 2 (m)s 2 1 , and its projection to Gr 2 ^ is 



A = 



~P3 ~ m 
P2 



1 
1 
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The r-function is 

t a = -(p2P3-Ea,2 +P3E1A + mE 2 ,4: +^£'2,3 - #3,4), 
where E^j :— (kj — k{) exp(#i + 9j). The contour plots of the solutions with m = and ra/0 are the 
same (except for some exponentially small difference) when i <C 0. In both cases, the plot consists of 
two line-solitons forming an X-crossing, where the parts of those solitons adjacent to the region with 
dominant exponential £^4 (i.e. for x 3> 0) are singular, see the left of Figure 10. 

On the other hand, for t 0, the contour plot with to = is topologically the same as it was for 
t -C 0, while the contour plot with to ^ has a box with dominant exponential £2,4, surrounded by 
four bounded solitons (some of which are singular). See the middle and right of Figure 10. So not only 
the contour plots but also the soliton graphs are different for t ^S> 0! 




Figure 10. The non-uniqueness of the evolution of the contour plots (and soliton 
graphs). The left panel shows the contour plot at t = —20 for any value of to. The 
middle panel shows the graph at t = 20 with to = 0, and the right one shows the graph 
at t — 20 with m = 1. These contour plots were made using the choice pi = 1 for all i, 
and (ki, . . . , K4) = (—2, —1, 0, 1.5). In all of them, the region at x has a positive 
sign (A3. 4 = 1) and other regions have negative signs. This means that the solitons 
adjacent to the region for i»0 are singular. 

Note that the non-uniqueness of the evolution of the contour plot (a tropical approximation) does 
not imply the non-uniqueness of the evolution of the solution of the KP equation as t changes. If one 
makes two different choices for the TOj's, the corresponding r-functions are different, but there is only 
an exponentially small difference in the corresponding contour plots (hence the topology of the contour 
plots is identical). 
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